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Abstract. We revisit a recent claim that the Earth's climate 
system is characterized by sensitive dependence to param- 
eters; in particular, that the system exhibits an asymmetric, 
large-amplitude response to normally distributed feedback 
forcing. Such a response would imply irreducible uncer- 
tainty in climate change predictions and thus have notable 
implications for climate science and climate-related policy 
making. We show that equilibrium climate sensitivity in 
all generality does not support such an intrinsic indetermi- 
nacy; the latter appears only in essentially linear systems. 
The main flaw in the analysis that led to this claim is in- 
appropriate linearization of an intrinsically nonlinear model; 
there is no room for physical interpretations or policy con- 
clusions based on this mathematical error. Sensitive depen- 
dence nonetheless does exist in the climate system, as well 
as in climate models — albeit in a very different sense from 
the one claimed in the linear work under scrutiny — and we 
illustrate it using a classical energy balance model (EBM) 
with nonlinear feedbacks. EBMs exhibit two saddle-node bi- 
furcations, more recently called "tipping points," which give 
rise to three distinct steady-state climates, two of which are 
stable. Such bistable behavior is, furthermore, supported by 
results from more realistic, nonequilibrium climate models. 
In a truly nonlinear setting, indeterminacy in the size of the 
response is observed only in the vicinity of tipping points. 
We show, in fact, that small disturbances cannot result in 
a large-amplitude response, unless the system is at or near 
such a point. We discuss briefly how the distance to the bi- 
furcation may be related to the strength of Earth's ice-albedo 



feedback. 
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1 Introduction and motivation 

1 . 1 Climate sensitivity and its implications 

Systems with feedbacks are an efficient mathematical tool 
for modeling a wide range of natural phenomena; Earth's 
climate is one of the most prominent examples. Stability 
and sensitivity of feedback models is, acc ordingly, a tradi- 



tional topic of theoretical c l imate studi es (Cess, 



197> 



6; 



Crafoord and Kalle 



Cess et al. 



1989) 



n, 



1978 



Schlesinger 



Roe an d Baker (2007) (RB07 hereafter) 



1985 



1976; Ghil, 



1986; 



have recently advocated existence of intrinsically large sen- 
sitivities in an equilibrium model with multiple feedbacks. 
Specifically, they argued that a small, normally distributed 
feedback may lead to large-magnitude, asymmetrically dis- 
tributed values of the system's response. 

Such a property, if valid, would have seriou s implications 
for climate dynamics dAllen and Frame , 



2007) and for mod- 
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eling of complex systems in general dWafkins and Freeman , 
2008). In this paper, we revisit the dynamical behavior of 
a general, equilibrium climate model with genuinely nonlin- 
ear feedbacks, and focus subsequently on a simple energy- 
balance model (EBM). We notice that the main, technical 
part of RB07's argument is well-k nown in the climate litera- 
ture, cf. ISchlesingen ( 1 19851 [1986), and thus it seems useful 
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to review the associated assumptions and possible interpreta- 
tions of this result. 

We rederive below in Section 11.21 the key equation of 
RB07 and comment on their purportedly nonlinear analysis 
in Section 11.31 We then proceed in Section [2] with a more 
self-consistent version of sensitivity analysis for a nonlin- 
ear model. This analysis is applied in Section [3] to a zero- 
dimensional EBM. Concluding remarks follow in Section|4] 

1 .2 Roe and Baker's (2007) linear analysis 

We follow here RB07 and assume the following general 
setup. Let the net radiation R at the top of the atmosphere 
be related to the corresponding average temperature T at the 
Earth's surface by R — R{T). Assume, furthermore, that 
there exists a feedback a = a(T), which is affected by the 
temperature change and which can, in turn, affect the radia- 
tive balance. Hence, one can write R = R (T, a(T)). 

To study how a small change A R in the radiation is related 
to the corresponding temperature c hange A T, one can use 
the Taylor expansion (lArfkenl 1 19851) to obtain, as A T tends 
to zero, 

Ai?=§|AT+||||AT + 0((AT) 2 ). (1) 

Here, 0(x) is a function such that 0(x) < C x as soon as 
< x < e for some positive constants C and e. 
Introducing the notations 

1 dR , , , dRda 
A a 1 a a a 1 

for the "reference sensitivity" Ao and the "feedback factor" 
/, we obtain 



AR= AT + ((AT) 5 



which readily leads to 



A 



AT = y— -Ai? + o ((AT) 2 ) , 



/ 

as long as / ^ 1. 

RB07 drop the higher-order terms in (O to obtain 

AT = j^-jAR, 



(2) 



(3) 



(4) 



which is their equation (S4). This equation leads directly 
to their main conclusion, namely that a normally distributed 
feedback factor / results in an asymmetric system response 



A T to a fixed forcing A R. The purported sensitivity is due 
to the divergence of the right-hand side (rhs) of Eq. [their 
equation (S4)] as / approaches unity. Figure[T] which is anal- 
ogous to Fig. 1 of RB07, illustrates this effect. 

Roughly speaking, RB07 use the following argument: If 
the derivative of R(T) with respect to T is close to 0, then 
the derivative of T with respect to R is very large, and a small 
change in the radiation R corresponds to a large change in the 
temperature T. Such an argument, though, is only valid for 
an essentially linear dependence R oc T. Our straightfor- 
ward analysis in Section [2] below shows that the sensitivity 
effect of Fig. 1 is absent in climate models in which gen- 
uinely nonlinear feedbacks R = R(T, a(T)) are present. 

It is worth noticing that, since one seeks the temperature 
change A T that results from a change A R in the forc- 
ing, it might be preferable to consider the inverse function 
T = T(R) or, more precisely, T = T(R, a(R)) and the 
corresponding Taylor expansion 

The main conclusions of our analysis will not be affected by 
the particular choice of direct or inverse expansion, provided 
the nonlinearities are correctly taken into account. 

1.3 Roe and Baker's "nonlinear" analysis 

In this section we address the analysis carried out by RB07 
in the supplemental on-line materials, pp. 4-5, Section "Non- 
linear feedbacks." The main conclusion of that analysis was 
that, given realistic parameter values for the climate system, 
the effects of possible nonlinearities in the behavior of the 
function R = R(T) are negligible and do not affect the sys- 
tem's sensitivity. We point out here two serious flaws in their 
mathematical reasoning that, each separately and the two to- 
gether, invalidate such a conclusion. 

First, and most importantly, despite their section's title, the 
analysis carried out by RB07 is still linear. Indeed, the Tay- 
lor expansion in their Eq. (S7) is given by 



ARk R'AT+ ^R" AT 2 , 



(5) 



where (■)' stands for differentiation with respect to T. But 
RB07 immediately invert this equation for A T subject to the 
assumption 

AT 2 = ATAT Q , 

where A To is a constant. Hence, instead of solving the 
quadratic Eq. |5j, Roe and Baker solve the following linear 
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approximation: 



AR 



(R' + -R" A T ) A T, 



and thus obtain the key formula [their Eq. (S8)] 

-A AR 



AT: 



l-/-|Ao R"AT ' 



(6) 



the last step uses the, correct, fact that R' = (1 — /)/Ao- 
This equation artificially introduces a divergence point for 
the temperature at / = 1 - (A /2) R"AT , which clearly 
cannot exist in a quadratic equation. Equation (O is thus a 
very crude approximation that significantly deviates from the 
true solution to the full quadratic equation (O — which we 
discuss below in Section |2] — and thus cannot be used to 
justify general statement s about climate models . 

The second flaw in the lRoe and Bakerl d2007l) reasoning is 



that, using the model and parameter values they suggest, one 
readily finds that: 

- R' = —2, i.e., global temperature and radiation are neg- 
atively correlated, which is hardly t he cas e for the cur- 
rent climate [e.g.. Held and Sodenl d2000l) l. We notice 
that the negative sign of the correlation follows directly 
from their Eq. (S10) and is not affected by particular 
values of the model parameters. Furthermore, 

- R" = —0.03, which means that the model they consider 
is, indeed, essentially linear, and thus not very realistic. 

Although, in this part of their analysis, Roe and Baker as- 
sumed that / = 0.4, it is easy to check that, for all / < 0.95, 
their model satisfies \R" /R'\ < 0.1 and is therewith very 
close to being linear. To conclude, the effects on nonlinear- 
ities are indeed negligible in the particular model studied in 
this part of the RB07 paper, since the model is very close to 
being linear; one cannot extrapolate, therefore, their conclu- 
sions to climate models with significant nonlinearities. 

We next proceed with a mathematically correct sensitivity 
analysis of a general climate model in the presence of truly 
nonlinear feedbacks. 



2 A self-consistent sensitivity analysis 

It is easily seen from the discussion in Section [L2l especially 
from Eq. (0, that the relationship (01 is a crude approxima- 
tion: it is valid only subject to the assumptions that (a) the 



higher-order terms in the expansion of A T are vanishingly 
small: 

0((AT) 2 ) < i— ^AT- 

and (b) the quantity in the rhs of this inequality is itself 
nonzero. 

If one assumes, for instance, that O ((AT) 2 ) ~ 
C(AT) 2 , where the precise meaning of g[x) ~ f(x) is 
given by lim x _,.o g(%) / f '{%) — 1> then the assumptions 
(a,b) above hold for A T that satisfy both of the following 
conditions 

0< AT« (l-/)/(A C) and 0<AT<e, (7) 

where C and e are defined after Eq. (HJ. The first of these 
conditions implies that the range of temperatures within 
which the approximation works vanishes as the feedback 
factor / approaches unity. Hence, all the results based on 
this approximation — including precisely the main conclu- 
sions of RB07 — no longer apply outside a vanishingly small 
neighborhood of / = 1.0. 

The asymptotic behavior we assumed above for 
O ((AT) 2 ) is not exotic. Consider for instance the 
function R = T 2 in the neighborhood of R = 0. Its Taylor 
expansion 

AR= 2TAT + ((AT) 2 ) 

can be used to obtain, ignoring the second-order term, 

AT ^ AR/(2T). (8) 

The last equation would seem to imply that the growth of 
A T is inversely proportional to T itself, so the change in T 
should increase infinitely fast as T goes to 0, a rather annoy- 
ing contradiction. 

The way out of this conundrum is to realize that the change 
A T given by Eq. (O is only valid in a small vicinity of T = 
and cannot be extrapolated to larger values. Of course, we all 
know that the function R = T 2 is nicely bounded and smooth 
in the vicinity of 0, but it is essential to take into account 
the second term in its Taylor expansion in this vicinity to 
obtain correct results. We show in Section [3]below that this 
simple example depicts the essential dependence of the Earth 
surface temperature on the global solar radiative input, for 
conditions close to those of the current Earth system. 

In summary, the linear approximation of the function 
R(T) derived by RB07 from its Taylor expansion is not valid 



4 



I. Zaliapin and M. Ghil: Climate Sensitivity 



when / approaches unity. In this case — which is precisely 
the situation emphasized by these authors — the higher-order 
terms "hidden" inside O ((AT) 2 ), which they neglected, are 
indispensable for a correct, self-consistent climate sensitivity 
analysis. 

A correct analysis of the case when / approaches unity 
needs to start with a Taylor expansion that keeps the second- 
order term 

AR = ^— -AT + a(AT) 2 + O ((AT) 3 ) , 

where a = R" /2. If O ((AT) 3 ) is much smaller than the 
other two terms on the rhs, then the temperature change can 
be approximated by a solution of the quadratic equation 



An 



-AT + a(AT) = AR. 



(9) 



The real-valued solutions to the latter equation, if they exist, 
are given by 



2 1 a An 



1-/ 

. a An 



4AR 



In particular, when / is close to 1 .0, then 

ATi, a «± 



AR 



(10) 



One can see from Eq. (TTOb that the proximity of the feed- 
back factor / to unity no longer plays an important role in 
the qualitative behavior of the equilibrium temperature. This 
point is further illustrated in Fig.|2]that shows the climate sys- 
tem's response AT as a function of the feedback factor / for 
different values of the nonlinearity parameter a. The most 
important observation is that the climate response does not 
diverge at / = 1; moreover, the asymmetry of the response 
due to the changes in feedback factor / rapidly vanishes as 
soon as the dependence of A R on A T becomes nonlinear. 

In general, one can consider an arbitrary number of terms 
in the Taylor expansion of R(T). The very fact that one re- 
lies on the validity of the Taylor expansion implies that R(T) 
is bounded and sufficiently smooth; in other words, a diver- 
gence of the equilibrium temperature due to a small change in 
the forcing contradicts the very assumptions on which RB07 
based their sensitivity analysis. 

3 Sensitivity for energy balance models (EBMs) 

We consider here a classical climate model with nonlinear 
feedbacks to illustrate that, in such a model: (i) the type of 



sensitivity claimed by RB07 does not exist; and (ii) sensitive 
dependence may exist, in a very different sense, namely in 
the neighborhood of bifurcation points, as explained below. 

3.1 Model formulation 

We consider here a highly idealized type of model that con- 
nects the Earth's temperature field to the solar radiative flux. 
The key idea on which these m odels are built is due to 
Budvkol (11969) and Sellers! (119691) . They have been subse- 
quently generalized and use d for many stud i es of c limate 



North 



stability a nd sensitivity ; see iHeld and Suarezl (119741) : 
J 19751) and bhill(ll976l) . among others. 

The interest and usefulness of these "toy" models resides 
in two complementary features: (i) their simplicity, which 
allows a complete and thorough understanding of the key 
mechanisms involved; and (ii) the fact that their conclusions 
have been extensively confirmed by studies using much more 
detailed and presumably realistic models, including general 
circulation models (G CMs); see, for instance, the reviews of 



North et al. 



(119811) and lGhil and Childress! (119871) . 
The main assumption of EBMs is that the rate of change of 
the global average temperature T is determined only by the 
net balance between the absorbed radiation R t and emitted 
radiation R a : 



of =A(T) 



Ro(T). 



(11) 



For simplicity, we f ollow here the zero - dimen sional 
(0-D) EBM version of Icrafoord and Kallenl d 19781) and 
Ghil and Childressl (11987b . in which only global, coordinate- 
independent quantities enter; thus 



Ri = nQo(l-a(T)), R = <rg(T)T I 



(12) 



In the present formulation, the planetary ice-albedo feed- 
back a decreases in an approximately linear fashion with T, 
within an intermediate range of temperatures, and is nearly 
constant for large and small T. Here Qq is the reference 
value of the global mean solar radiative input, a is the Stefan- 
Boltzmann constant, and g(T) is the grayness of the Earth, 
i.e. its deviation from black-body radiation aT i . The pa- 
rameter [i ss 1.0 multiplying Q indicates by how much the 
global insolation deviates from its reference value. 
We model the ice-albedo feedback by 

(t \ r i-tanh(K(r-r c )) 

a(T;K) — ci + c 2 . (13) 
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This parametrization represents a smooth interpolation be- 
tween the p iecewise-lin ear formula of Sellers-type mo dels, 
like those of bhill J 19761) or lcrafoord and Kallenl dl978l) . and 



the piec ewise-constant fo r mula of Budyko-type models, like 



Held and Suarezl (1 1974 or 



North 



(1975). 



those of 

Figure [3^ shows four profiles of our ice-albedo feedback 
a(T) = a(T) as a function of T, depending on the value 
of the steepness parameter k. The profile for k> 1 would 
correspond roughly to a Budyko-type model, in which the 
albedo a takes only two constant values, high and low, de- 
pending on whether T < T c or T > T c . The other profiles 
shown in the figure for smaller values of k, correspond to 
Sellers-type models, in which there exists a transition ramp 
between the high and low albedo values. Figure shows 
the corresponding shapes of the radiative input il; = i?;(T). 
The greenhouse effec t i s parametrized, as in 



Crafoord and Kallenl ( 1 19781) and lGhil and Childress! ( 1 1987b 
by letting 



g(T) = 1 - m tanh ((T/T ) 6 ) 



(14) 



Substituting this greenhouse effect parametrization and the 
one for the albedo into Eq. ( fTTb leads to the following EBM: 



cT = fi Q (l-a(T)) 

- aT i [1 - m tanh ((T/T ) 6 )] 



(15) 



where T — dT/dt denotes the derivative of global tempera- 
ture T with respect to time t. 

It is important to note that current concern, both 
scientific and public, is mostly with the greenhouse 
effect, rather than with actual changes in insolation. 
But in a simp le EBM model — whether globally aver- 
aged, like in ICrafoord and Kallen (11978b and here, o r 



coordinate-dependent, a s in iBudvkol dl969b":ISellersl d!969b ; 



Held and Suarezl dl974h ; iNorthl d 1975b or |Ghil| d 1976b — 
increasing /.i always results in a net increase in the radiation 
balance. It is thus convenient, and quite sufficient for the 
purpose at hand, to vary fi in the incoming radiation R[, 
rather than some other parameter in the outgoing radiation 
R D . We shall return to this point in Section[4] 

3.2 Model parameters 

The value S of the solar constant, which is the value of 
the solar flux normally incident at the top of the atmosphere 
along a straight line connecting the Earth and the Sun, is as- 
sumed here to be S — 1370 Wm _1 . The reference value of 



the global mean solar radiative input is Qq = 5/4 = 342.5, 
with the factor 1/4 due to Earth's sphericity. 

The parameterization of the ice-albedo feedback in 
Eq. O assumes T c = 273 K and ci = 0.15, c 2 = 0.7, 
which ens ures t hat a (T) is bounded between 0.15 and 



0.85, as in Ghi]| dl976b : see Fig. [3^. The greenhouse effect 
parametrization in Eq. dl~4l > uses m = 0.4, which corre- 
sponds to_40% clpudcover, and T ~ 6 = 1.9 X 10~ 15 K~ 6 



(Sellers, 



1969; 



Ghil, 



is a 



1976). The Stefan-Boltzmann constant 
5.6697 x 10~ 8 Wm~ 2 K~ 4 . 



3.3 Sensitivity and bifurcation analysis 

3.3.1 Two types of sensitivity analysis 

We distinguish here between two types of sensitivity analy- 
sis for the 0-D EBM (fTTb . In the first type, we assume that 
the system is driven out of an equilibrium state T = Tq, 
for which Ri(Ta) — R (Tq) = 0, by an external force, and 
want to see whether and how it will return to a new equi- 
librium state, which may be different from the original one. 
This analysis refers to the "fast" dynamics of the system, and 
assumes that Ri(T) - R D (T) / for T / T ; it is of- 
ten referred to as linear stability analysis, since it consid- 
ers mainly small displacements from equilibrium at t = 0, 
T(Q) = T + 6(0), where 0(0) is of order e, with < e < 1, 
as defined in Section[2] 

The second type of analysis refers to the system's "slow" 
dynamics. We are interested in how the system evolves 
along a branch of equilibrium solutions as the external force 
changes sufficiently slowly for the system to track an equi- 
librium state; hence, this second type of analysis always as- 
sumes that the solution is in equilibrium with the forcing: 
Ri(T) - R (T) = for all T of interest. Typically, we want 
to know how sensitive model solutions are to such a slow 
change in a given parameter, and so this type of analysis is 
called sensitivity analysis. In the problem at hand, we will 
study — again followin 



Ghil and Childress 



ing 

37) 



Crafoord and Kallen 



dl978b and 



1987) — how changes in /i, and hence 
in the global insolation, affect the model's equilibria. 

A remarkable property of the EBM governed by Eq. (fTTT > 
is the existence of several station ary solutions that describ e 



equilibrium climates of the Earth ( Ghil and Child ress. 



1987) 



The existence and linear stability of these solutions result 
from a straightforward bifurcation analysis of the 0-D EBM 
dTTJ, as well as of its one-dimensional, latutude-dependent 
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counterparts (Ghil, 1976. Il994l) : there are two linearly sta- 
ble solutions — one that corresponds to the present cli- 
mate an d one that correspond s to a much colder, "snowball 
Earth" ([Hoffman et all [1998) — separated by an unstable 
one, which lies about 10 K below the present climate. 

The existence of the three equilibria — two stable and 
one unstable — has been confirmed by such results be- 
ing obtained by several dist i nct EBMs, of eithe r Budyko- 

Nonlin- 



Sellers-type jNorth et all ll98lL bhili [l99J). 



ear stability, to large perturbations in the initial state, has 
been investigated by introducing a varia tiona l princ iple for 



1976) and of 



the latitu de-dependen t EBMs of Sellers (Ghil, 
Budyko (INorth et all 1 1 98 If) type, and it confirms the linear 
stability results. 

3.3.2 Sensitivity analysis for a 0-D EBM 

We analyze here the stability of the "slow," quasi-adiabatic 
(in the statistical-physics sense) dynamics of model (fl31 l. 
The energy -balance condition for steady-state solutions R[ = 
R takes the form 

MQo (1 - a(T)) = (tT 4 [1 - m tanh ((T/T ) 6 )] . (16) 

We assume here, following the previously cited EBM work, 
that the main bifurcation pa rameter is ii\ this hap pens to 
agree with the emphasis of iRoe and Baker! (120071) on cli- 
mate sensitivity as the dependence of mean temperature T 
on global solar radiative input, denoted here by Q = fiQo- 

Figure [4] shows the absorbed and emitted radiative fluxes, 
i?i and R , as functions of temperature T for /i = 0.5, 1 and 
2.0. One can see that Eq. ([TBI may have one or three solu- 
tions depending on the value of /i: only the present, relatively 
warm climate for /i = 2.0, only the "deep-freeze" climate for 
/U = 0.5, and all three, including the intermediate, unstable 
one for present-day insolation values, /i = 1.0. These steady- 
state climate values are shown as a function of the insolation 
parameter /.i in the bifurcation diagram of Fig. [5] 

The "fast" stability analysis (not presented here) shows 
that small deviations 6(0) from an equilibrium solution, 
while all parameter values are kept fixed, may result in 
two types of dynamics, depending on the initial equilibrium 
To: fast increase or fast decrease of the initial deviation 



dGhil and Childressl 1 19871) . The fast increase characterizes 
unstable equilibria: a small deviation 8(0) from such an equi- 
librium To forces the solution to go further and further away 
from the equilibrium. In practice, such equilibria will not be 



observed, since there are always small, random perturbations 
of the climate present in the system: just think of weather as 
representing such perturbations. 

The fast decrease of the initial deviation 9(0) character- 
izes stable solutions; only such equilibria can be observed in 
practice. The two stable solution branches of (l5[ are shown 
by solid lines in Fig. [5] while the unstable branch is shown by 
the dashed line. The arrows show the direction in which the 
temperature will change when drawn away from an equilib- 
rium by external forces. This change, whether away from or 
towards the nearest equilibrium, is fast compared to the o ne 
that occurs along either solution branch ( IGhill 1 1 97611 1 9941) . 



3.3.3 Bifurcation analysis 

Given the choice of model parameters, the present climate 
state corresponds to the upper stable solution of Eq. ( fT5l ). 
at fi = 1 (see Fig. [5]). It lies quite close to the bifurcation 
point (/i, T) « (0.9, 280 K), where the stable and unstable 
solutions merge. 

The so-called normal form of this bifurcation is given by 
the equation 

X = fl-X 2 , (17) 

where X is a suitably normalized form of T, and Jx is a nor- 
malized form of fj,. Equation ( TPTT i describes the dependence 
between T and \i in a small neighborhood of the bifurca- 
tion point. In particular, the stable equilibrium branch is de- 
scribed by 

X = 

this result has exactly the same form as the positive solution 
of Eq. ( [Tol l, given by our self-consistent analysis of climate 
sensitivity in the presence of genuine nonlinearities, cf. Sec- 
tionlU Hence, the derivative dX/dp,, and thus dT/d/i, goes 
to infinity as the model approaches the bifurcation point; this 
is exactly the situation discussed earlier in Section|2] 

It is important to realize that the parabolic form of 
temperature dependence on insolation change is not an 
accident due to the particul arly simple form of EBMs. 
Wetherald and Manabd (11975b clearly showed, in a slightly 
simplified GCM, that not only the mass-weighted temper- 
ature of their total atmosphere, but also the area-weighted 
temperatures of each of their five model levels, exhibits such 
a parabolic dependence on fractional radiative input; see 
Fig. 5 in their paper. Moreover, these authors emphasize that 
"As stated in the Introduction, it is not, however, reasonable 
to conclude that the present results are more reliable than 



I. Zaliapin and M. Ghil: Climate Sensitivity 



7 



the results from the one-dimensional studies mentioned 
above simply because our model treats the effect of transport 
explicitly rather than by parameterization. [...] Nevertheless, 
it seems to be significant that both the one-dimensional and 
three-dimensional models yield qualitatively similar results 
in many respects." 

In fact, rigorous mathematical results demonstrate that 
the saddle-node bifurcation whose normal form is given 
by Eq. ( ITTb occurs in several systems of nonlinear partial 
differentia l equat i ons, such as the Na vier-Stokes equations 



dConstantin et al. 



1989; 



Temam, 



1993), and not only in ordi- 



nary differential equations, like Eqs. dTTb and (fT5T l above. We 
emphasize, though, that this does not cause the temperature 
to increase rapidly due to small changes in insolation: the 
presence of the bifurcation point will result in small, positive 
changes of global temperature for slow, positive changes in 
[i, while it may throw the climate system into the deep-freeze 
state for slow, negative changes in fi. 

4 Discussion 

4.1 How sensitive is climate? 

Making projections of climate change for the next decades 
and centuries, evaluating the human influence on future Earth 
temperatures, and making normative decisions about cur- 
rent and future anthropogenic impacts on climate are enor- 
mous tasks that require solid scientific expertise, as well as 
responsible moral reasoning. Well-founded approaches to 
handle the moral aspects of the p roblem are still being de- 



bated [e.g., iHillerbrand and Ghill ( 120081) 1. It is that much 



more important to master existing tools for acquiring accu- 
rate and reliable scientific evidence from the available data 
and models. Several of these tools come from the re alm of 

)renzl \l9< 



Smalel 


196' 


• 


Ruelle and Takens 




1971; 


1987 


Ghil, 


1994; 


Ghil et al. 


2008 


). 



Ghil and Childress, 



A straightforward analysis, carried out in Section [2] of 
this paper, shows that a proper treatment of the higher-order 
terms in a climate model with nonlinear feedbacks does not 
reveal the exaggerated sensitivity to forcing that was used in 
RB07 to advocate intrinsic unpredictability of climate pro- 
jections. We emphasize that the error in Roe and Baker's 
analysis is not related to their choice of model formulation or 
of the model parameters nor to their interpretation of model 
results. The problem is purely a matter of elementary calcu- 



lus, and is due to inappropriate, and unnecessary, lineariza- 
tion of a nonlinear model. 

Our analysis complemen ts, reinforces and goes beyond 



that of 



Hannart et al. 



( 2009), who also showed that the claim 
of RB07 "results from a mathematical artifact." We notice 
simply that Hannart and colleagues did not even question the 
linear approximation framework of RB07 and still concluded 
that the claims of irreducibility of the spread in the envelope 
of climate sensitivity are not supported by the RB07 analysis. 

To summarize, while the general human concern about cli- 
mate sensitivity expressed by RB07 should be reasonably 
shared by many, their scientific conclusions do not follow 
from their model and its results, when correctly analyzed, as 
done here in Section|2] Nor are these conclusions supported 
by other models of greater detail and realism, when properly 
investigated. Accordingly, conclusions about the likelihood 
of extreme warming resulting from small changes in anthro- 
pogenic forcin g can hardly be used t o support political pro- 
posals \e.s.. lAllen and Frame! ( 120070 1 that claim to provide 
future directions for the climate-related sciences. 

Still, this paper's analysis does not preclude in any sense 
the Earth's temperature from rising significantly in com- 
ing years. The methods illustrated here can only be used 
to study climate sensitivity in the vicinity of a given state; 
they cannot be applied to investigate climate evolution over 
tens of years, for example in response to large increases 
in greenhouse gases or to other major changes in the forc- 
ing, whether natural or anthropogenic. This latter prob- 
lem requires global interdisciplinary efforts and, in partic- 
ular, the analysis of the entire h ierarchy of climate models 
dSchneider and Dickinson! 1 19741). from conceptual to inter - 



mediate to fully coupled GCMs (Ghil and Robertson, 



2000) 



It also requires a much more c areful study of ran dom effects 



than has been done heretofore (Gh il et al. 



2008). 



It seems to us that Roe and Baker's title question "Why Is 
Climate Sensitivity So Unpredictable?" still remains open. 

4.2 Where are the "tipping points"? 

The S-shaped diagram o f Fig. |5] — s ee also Fig 10.6 in 



Ghil 



(1994) — was 



Ghil and Childress! 11987) and Fig. 4 in|| 
used here to show the smoothness and boundedness of tem- 
perature changes as a function of insolation changes, away 
from a saddle-node bifurcation, like that of Eq. (TTOb in Sec- 
tion|2]or of Eq. (QjJ in Section [3331 
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This S-shaped curve nevertheless reveals the existence of 
sensitive dependence of Earth's temperature on insolation 
changes, or on other changes in Earth's net radiation budget, 
such as may be caused by increasing levels of greenhouse 
gases, on the one hand, or of aerosols, on the other. This sen- 
sitive dependence is quite different from the one advocated 
by RB07. Namely, if the parameter [i were to slightly de- 
crease — rather than increase, as it seems to have done since 
the mid-1970s, in the sense described in the last paragraph 
of Section [3TTI — then the climate system would be pushed 
past the bifurcation point at /i w 0.9. The only way for the 
global temperature to go would be down, all the way to a 
deep-freeze Earth, with much lower temperatures than those 
of recent, Quaternary ice ages. 

It has become common in recent discourse about po- 
tentially irre versible clima t e cha nge to talk about "tipping 

(12008b . The term was o riginally 



Lenton et al. 



points ; e.g. 

introduced into the social sciences by lGladw ell (2000) to de- 
note a point at which a previously rare phenomenon becomes 
dramatically more common. In the physical sciences, it has 
been identified with a shift from one stable equilibrium to an- 
other one, i.e., to a saddle-node bifurcation, as seen in Fig. [5] 
here and explained in Section [3.3.3l above. 

In the EBM context of Fig. [5] it would require an enor- 
mous, almost twofold increase in the insolation in order for 
a deep-freeze-type equilibrium to reach the bifurcation point 
at fi rj 1.85 and jump from there to T f=s 350 K, a tem- 
perature that sounds equally unpleasant. Within the broader 
context of the recent debates on how to exit a snowball-Earth 
state, very large, and pos sibly implausible incr eases in CO2 
levels would be required dPierrehumbertt 120041) . 

Indeed, the likelihood to actually reach the tipping point 
to the left of the current climate in Fig. [5] seems to be quite 
small. Mechanisms for entering a snowball-Earth climate 
have been recen tly studied with a number of fairly r ealistic 



Poulsen and Jacob 



climate models (Hyde et al. , 2000; Donnadieu et al. 



2004 



2004). Both modeling and independent 
geological evidence suggest that Earth's climate can sustain 
significant fluctuations of the solar radiative input, and hence 
of global temperature, without entering the snowball Earth, 
and evidence for Earth ever having been in such a state is still 
controversial. 

Nevertheless, the existence of the upper-left tipping point 
shown in Fig. [5] is confirmed by numerous model studies, 
including GCMs, and we have already cited some evidence 
also for the lower-right tipping point in the figure. Several 



hypothetical tippin g points on the "warm" side have been 
identified by 



Lenton et al. 



(2008) and references therein, 



among many others. But only few of these have been stud- 
ied with the same degree of mathematical and physical de- 
tail as the ones of Fig. [5] here. One worthwhile example is 



lation ( 


Stommel, 


1961 


Brvan, 


1986; 


Ouon and Ghil, 


1992; 


Thual and McWilliams 


1992:1 


Diikstra and Ghil, 


2005). 



5UTSU- 



2006; 



Accordingly, humankind must be caref ul — in 
ing its recent inter est in geoengineering dCrutzen , 
MacCrackenl |2006) — to stay a course that runs between 
tipping points on the warm, as well as on the "cold" side of 
our current climate. In any case, the existence, position and 
properties of such tipping points need to be established by 
physically careful and mathematically rigorous studies. The 
"margins of maneuver" seem reasonably wide, at least on 
the time scale of tens to hundreds of years, but this does not 
eliminate the possibility to eventually reach one such tipping 
point, and thus we are led directly to the next question. 

4.3 How close are we to a cold tipping point? 

Let us assume for the moment that the dangers of further 
warming will lead humanity to actually stop, and possibly 
reverse, the current trend of an increasingly positive net ra- 
diation balance. Given, on the other hand, the dangers of a 
snowball Earth, one might want to estimate then the close- 
ness of the climate system to the top-left bifurcation point in 
Fig.[5]here. 

The GCM simulations of 



Wetherald and Manabe (1975) 



(see again their Fig. 5) suggest that this point might lie no 
farther than 5% below the current value of the solar constant. 
At the same time, the Sun has been much fainter 4 Gyr ago 
(by approximatyely 25-30%) than today, without the Earth 
ending up in a deep freeze, except possibly much later. So 
how close are we to this tipping point? 

Figure [6] here shows stable and unstable equilibrium solu- 
tions for different profiles of the ice-albedo feedback, a = 
a(T; k); this profile is determined by the value of the steep- 
ness parameter k (cf. Fig. [3^). The figure suggests that the 
steeper the ramp of the ice-albedo feedback function, i.e. the 
larger k, the further away the bifurcation might lie. In fact, 
for a very smooth dependence of the albedo on temperature, 
i.e. for a very small k, there is no bifurcation at all (not 
shown): very small values of k produce only a single-valued, 
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smoothly increasing, stable equilibrium solution to $15[ for 
any value of /i. 

It seems worthwhile to carry out systematic bifurcation 
studies with atmospheric, oceanic and coupled GCMs to ex- 
amine this question more carefully, for warm tipping points, 
as well as for cold ones. Such studies are made possible by 
current computing capabilities, along wi th well-developed 
methods of numerical bifurca tion theory (IDiikstra and Ghil , 



2005 



Simonnet et al. 



2009). This approach holds some 



promise in evaluating the distance of the current climate state 
from either a catastrophic warming or a catastrophic cooling. 
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Fig. 1. Temperature response AT as a function of the feedback 
factor / in a linear climate model; modified after iRoe and Bakei 
1 2007). In the absence of feedbacks (/ = 0), the response is given 



by A T — A To = 1.2°C. The response is amplified by feedbacks; 
it diverges (AT — 5* oo) as / —5- 1. The shaded areas illustrate the 
hypothetical symmetric distribution of the feedback factor / and the 
corresponding asymmetric distribution of the system response A T. 
Our study shows that this prominently asymmetric response is only 
seen in a linear model and is absent in a general nonlinear model 
(see Fig.[2]below). 




0.2 



0.4 0.6 
Feedback factor, f 



0.8 



Fig. 2. Temperature response A T as a function of feedback factor 
/ in a nonlinear, quadratic climate model, governed by Eq. l[9}. The 
curves correspond (from top to bottom) to increasing values of the 
magnitude a = R" /2 of the quadratic term; the values of a for each 
curve are given in the figure's legend, in the upper-left corner. The 
uppermost curve (red) corresponds to a = and is the same as the 
one shown in Fig. Q] above. Extreme sensitivity, expressed in the 
divergence of A T, is only seen in the linear model; it rapidly van- 
ishes in the nonlinear model, as the nonlinearity factor a increases. 
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Fig. 3. Dependence of the absorbed incoming radiation Ri on the 
steepness parameter k: (a) ice-albedo feedback a = a(T; k), and 
(b) absorbed radiation Ri = Ri(T; k), for different values of k; see 
Eqs. 03 and (B). 
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Fig. 4. Outgoing radiation R a (red, dashed line) and absorbed in- 
coming solar radiation Ri (blue, solid lines) for our 0-D energy- 
balance model (EBM), governed by Eq. H5\ . The absorbed radia- 
tion is shown for n — 0.5, 1 and 2.0 (from bottom to top), while 
k = 0.05. Notice the existence of one or three intersection points 
between the _R curve and one of the Ri curves, depending on the 
value of these points correspond to the equilibrium solutions of 
dl5t . i.e. to steady-state climates. 
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Fig. 5. Equilibrium solutions of the EBM d 1 5b depending on the 
fractional change /j, in insolation. Notice the existence of two sta- 
ble (blue, solid lines) and one unstable (red, dashed line) solution 
branches. The arrows show the direction in which the global tem- 
perature will change after being displaced from a nearby equilib- 
rium state by external forces. The current Earth state corresponds 
to the upper stable solution at /i = 1; in this figure k = 0.05. 
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Fig. 6. Equilibrium solutions of the EBM J 1 5 b for ice-albedo feed- 
back functions a — a(T; k) given by Eq. J 1 3 1 > and corresponding 
to the steepness parameter values k = 0.01, 0.05, 0.1 and 1. Notice 
the deformation of the stable (blue, solid lines) and unstable (red, 
dashed lines) solution branches: the "cold" tipping point moves to 
the left, away from fj, — 1.0, as k increases; compare Fig.[3ji. 



